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Abstract: 

The universality of correlation functions of eigenvalues of large random matrices has 
been observed in various physical systems, and proved in some particular the 
hermitian one-matrix model with polynomial potential. Here, we consider the more difficult 
case of a unidimensional chain of matrices with first neighbour couplings and polynomial 
potentials. 

An asymptotic expression of the orthogonal polynomials allows to find new results for 
the correlations of eigenvalues of different matrices of the chain. 

Eventually, we consider the limit of the infinite chain of matrices, which can be in- 
terpreted as a time dependent one-matrix model, and give the correlation functions of 
eigenvalues at different times. 
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It has now since a long time been observed experimentally or numerically |l],0,|3| that 
the distribution of energy levels of disordered systems is universal in some regime. For 
instance the connected correlation function of two levels separated by a small number of 
other levels does not depend on the details of the system, only on its symmetries, while 
the density of levels is very dependent of the specific details of the system. It was thus 
conjectured that the correlation functions can be obtained from a gaussian random matrix 
model (the matrix might be the hamiltonian, the scattering or transfert matrix). Such a 
conjecture would be a kind of a central-limit-theorem for large random matrices. 

This conjecture has been proved in the special case of a one hermitian matrix model 
with polynomial potential and for a two-matrix- mo del ||. It has also been noted that 
the connected correlation functions of more than two eigenvalues should present a stronger 
universality than the density itself. 

The analysis of was based on the method of orthogonal polynomials. The cor- 
relation functions were expressed in terms of kernels, depending on two variables, which 
are given as sums of polynomials. Those results are exact and have been known for a long 
time. The problem was to find an asymptotic expansion in the large iV limit (N is the size 
of the matrices), indeed those kernels involve the sum of N polynomials of degree running 
from to N — 1. In the one- matrix model case [|J, the Darboux-Christoffel theorem al- 
lows to rewrite the kernel with only two polynomials of degree N and N — 1. Asymptotic 
expressions of orthogonal polynomials are then used to evaluate the universal correlation 
functions in the short range regime. 

In H it was claimed that this analysis could probably be extended to a more general 
case, which is the chain of p random hermitian matrices Mi, . . .M p , where each matrix 
Mk is coupled linearly to the following one Mk+i- 

In particular, when the number of matrices of the chain becomes infinite and the 
coupling are chosen appropriately, this model can be viewed as a time-dependent random 
matrix. The coupling between neighbouring matrices of the chain being then a kinetic 
term of the form (dM /dt ) 2 . 

Here, we will generalize the analysis of to the chain of matrices. The paper is 
organized as follows: 

The first section concerns the discrete chain, and the second section, the continuous- 
time limit. In the first section, we first present the matrix- model, remind the orthogonal 
polynomial method, then we relate the correlation functions to the orthogonal polynomials 
via the kernels and generalize the Darboux Christoffel theorem in order to rewrite those 
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kernels as a sum of a finite number of terms. A WKB approximation of the orthogonal 
polynomials allows to find asymptotic expressions of the kernels, and thus to find the 
correlation functions in the large iV limit. We then conclude by examining the universal 
properties of those correlations. 



1. The Chain of matrices 

Let us first present the model and introduce the notations coherent with those of 0. 
Consider a linear chain of p + 1 random hermitian N x N matrices Mj (0 < i < p) , 
with a probability: 

1 JL -NtrVdMi) -Nctr (M, - M i+l ) 2 , s 

V(M l ) = -l[e XT e ( L1 ) 

i=0 i=0 

where the Vi are polynomial potentials, c is the coupling constant between neighbouring 
matrices, and Z is the partition function. ( In the next section, we will consider the 
continuum limit of this model: the index i will become a continuous variable: the time 
t = ze, and with 2c = 1/e, the linear term £) f ctr (Mi — Mi+i) 2 will become a kinetic term 
JiM 2 (t)dt ). 

The Itzykson-Zuber formula allows to integrate out the angular variables (the 
unitary group), and leaves us with the joint probability for the eigenvalues ( let us note 
Xi-j (0 < j < N — 1) the j eigenvalue of the matrix Mi): 

piXi-j) = -A(Ai)A(A p ) J]e 11 detfc ^ e ( L2 ) 

i=0 i=0 V / 

where the A are the Vandermonde determinants: 

A(X l ) = l[(X i;k -X i;l ) 

k<l 

We would now like to compute the conditional probabilities of some subset / of these 
iVx (p + 1) eigenvalues. We thus have to integrate (|1.2| ) over all the eigenvalues which do 
not belong to /. For instance, the density of the eigenvalues of Mi is: 




the correlation function of two eigenvalues of Mj is: 

P«(Aj ; i, A i;2 ) = / p ][[ dA j;fc 
Cj,fc)/(«,i),(*.2) 

and the correlation function of two eigenvalues of two matrices Mj and Mj is: 



Pij(K;lAj;l) = J P II dA 



Z;fc 

(Lfc)/(i,l),(j,l) 

As in the one matrix case all these densities and correlation functions can be 

calculated by the orthogonal polynomials method ||, let us recall this method M. 

1.1. Orthogonal polynomials 

Consider two families of polynomials V n and V n , of degree n, beginning with the same 
leading term, and which obey the orthogonality relation: 



/ 



dA ...dA p e e V n {\o)V m {X P ) = 5 n , m (1.3) 



we define the wave functions ip n and by: 

^n(Ao) = V n (X )e 



-N\V P (\ P ) 

i>n(X P ) = V n (Xp) e 

(note that the normalizations differ from ||). with the help of the orthogonality relation 
( |1.3| ) , we can define two families of Hilbert Spaces Ei,Ei, and orthogonal functions in 
each of them: 



V'O.nCAo) =V , n(A ) 

-iVc((A -A 1 ) 2 + ... + (A l _ 1 -A l ) 2 ) 



i>i;n{Xi) = / dA . ..dAj_i ip n (\ )e 



-iV(|Vb(A ) + Vi(Ai) + . . . + Vi-^Xi-i) + iVi(Xi)) , s 
e - - (1.4) 



^p;n(A p ) =$ n (\ p ) 

7 (\^ [a\ hx 7 r\ \ -^((A i -A, +1 ) 2 + ... + (A p _ 1 -A p ) 2 ) 
Wi;n{^i) = / dA i+i . . . dA p ip n {\ p ) e 

-N^ViiXi) + V l+1 (X l+1 ) + ... + V p _i(Vi) + I^p(Ap)) , . 
e - " (1.5) 



We shall note them with the convenient Dirac notations: 

<n\i = ip i;n , \n>i=^i- n . (1.6) 

In the space Ei, we have the orthogonality relation: 

/dA j * ! „(A i )^ m (A i )=< n | m > i= S n ,„ 

In each of these spaces, we can define the usual operators (acting on the right hand side: 
the ket): 

• Aj, the operator which multiplies ipi(Xi) by Aj. 

• Pi = which differentiate ipi with respect to A,. 

These operators are defined only in the Hilbert space Ei. However all the Ei are 
isomorphics, and an operator O initially defined in Ei can also be defined in any Ej by its 
matrix elements: 

< n\6\m >j = f < n\6\m >i= J d\i ipi ;n (K) O ^ ;m (Aj) 

From now on, we will drop the index i for the bras and kets. 

1.2. Equations of motion 

From the former definitions we immediately obtain the equations of motion: 

Pi = 2c(X i+1 - Xi) - ^V-(Xi) (1.7) 

P = 2c(A 1 -Ao)-^V / (A ) , P p = 2cCX p -X p - 1 ) + ^V;(X P ) (1.8) 
and with an integration by parts: 

A i _ 1 + A i+1 -2A i = ^-V; , (A i ) (1.9) 

2c 

Let us now introduce more convenient notations. Since we began with polynomials V n 
and V n , we know how the multiplications or derivations by Ao or A p will act: multiplication 
by A p raises the degree of V n by 1, and X p V n (X p ) can be decomposed onto the base of the 
V n -k with k > —1: 

X p V n (X p ) = a(n)V n+1 (X p ) + ^2a k (n- k)V n - k (X p ) 

k>0 
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(where a(n) is the ratio of the leading coefficients of V n and V n -\-i, and the ak(n) are 
coefficients to be determined later). 

Let us write this in operatorial notations. In this purpose we introduce the operator 
x which decreases the level ( a kind of anihilation operator) : 

x\n >= \n — 1 > , < n\x =< n + 1| 

Although x is not invertible, we shall abusively write x^ = i, for it will make no difference 
when we go to the large N limit, and it will considerably simplify the notations *. 
We can then write: 



Aq =\a{h) + x k ak{h) 

k>0 

A p = — a(n) + ^ x k ak(n) 



x 

k>0 

(1.10) 



x k>0 

Remember that Ao acts on the ket \n >, i.e. on the polynomial V n , its adjoint acts on the 
bra V n . Note also that the first term, a(n) is the same for both A and A p because we 
have chosen the polynomials V n and V n with the same leading coefficient. 

Similarly, noting that the operator P = decreases the degree of the polynomial 

V n (X p ) by 1: 

-tt— V n (X p ) = — — nV n -i + ... 

dX p a(n — 1) 

we can can express the operators P$ and P p in power series of x: 

i i f> 

^ + 2W = ^iV +0(i2) (L12) 
We might as well write any of the operators A^ with such notations: 

K = ^2ai. k (h)x k 7 
k 

but let us first go to the large iV limit. 



* we have x|0 >= 0, and x 1 is not denned only on |0 >. One solution could be to define a 
state | — 1 >, provided that all the «fe(— 1) vanish, which is true. 
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1.3. large N limit 

In the classical limit N — > oo, all those operators become numbers. Indeed, the 
commutators [P, A] = 1/N and [x,n/N] = x/N are proportional to 1/N which thus plays 
the role of Ti **. We then write: 

K(x) = J2 a i,kX k with 11 degVj <k< n de S^i ( L13 ) 

k j=0 

The bounds on k are easily derived from the equations of motion and boundary conditions. 
We also consider the limit where n is large and close to N, so that the don't depend 
anymore on n, they are just coefficients. 

In addition, There exist a remarkable relation (the proof from the cannonical commu- 
tation relations is not difficult but of no interest for what follows) : 

1 = c^/c(a i+ i ifc - ai-i )k )ai-k- (1-14) 
k 

Let us rewrite in the classical limit the equations of motion (|1.9| ) and boundary conditions 
( p..ll|) (|1.12|) previously written for operators. We have the following system of equations: 

A l _ 1 +A l+1 -2A l = ^-y/(A l ) l<i<p-l (1.15) 

with the boundary conditions: 

y '(A ) - 2c(Ai - Ac) = — + 0(l/x 2 ) (1.16) 

ax 

V;{\) - 2c(A p _! - A p ) = ^ + 0(x 2 ) (1.17) 

One can verify that we have exactly as many equations as unknowns. If we were able to 
solve this system of algebraic equations and determine all the a^k, we could define p + 1 
functions \i(x), of an auxiliary variable x. We will see below the important role they play. 

** actually, this is true only if the support of the density is connected, i.e. we assume we have 
a one-cut solution, for a fc-cut solution, we would need to consider the operators as k x k matrices, 
e.g. for a symetric double well, one needs to distinguish between even and odd values of n, which 
introduces two sets of coefficients ak(2n) and otk(2n + 1). 
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1.4- WKB approximation 

One can find (by a simple generalization of |5J, i.e. by a kind of saddle point method 
for matrix integrals to the explicit expressions given in appendix B) some asymptotic 
expressions of the in the limit N large and \n — N\ ~ 0(1): 

i r\\ XT ( 71 ^" p/4 1 1 n—N -2Ncf(X 1 -X )y fVb(Ao) 

x I A (x)=Ao v uv ' 

, { \ \ ^/7T V /2-p/4 1 1 n _ N -Ncf'jXi+i-Xi-M 

, /-a \ WM p/4 1 1 n-N -2Ncf x (X p - Ap_i)Ai -f^ p (A p ) 
W A p ~ V tt -7= 7 g e e 

7 ,n WM p/4 1 1 N-n-i 2NcJ x (X 1 -X )X' -fV (X ) 



1.18) 



X 



2 " J-K( 



X 



We shall not prove those asymptotic expressions, but just give some intuitive explanations. 

- First, observe that at leading order, all of them have the form: 

-Nf Xi PidXi j ~ Nf Xi PidXi 

ipi ~ e and ipi ~ e 

which is simply the solution of the differential equation Pi\n >= jj-^r\n > in the large 
N limit . The lower bound of the integral, which has not been written here for simplicity, 
depends on i, it is determined by the condition that j Xl I^'(//j)d/ij = T^(Aj) . 

- The x n ~ N term comes from the definition of x: 

< N\x n ~ N =< n\ 

- Moreover, observe that the approximation for V'i+l can be derived from ^ by steepest 
descent in (1.4 ), and the expressions for the ifi's can be derived from the ifi's by x — > 1/x 
and i — ► p — i. 



- Finally, the normalization constants and the 1 / ^2'irX' i (x) are just what is needed to 
satisfy the normalization condition 

< n\m >= 5 nm . 

Remark that all that is nothing else than WKB approximation. 

Remember that, in quantum mechanics, the wave function of a particle outside a 
potential well decreases exponentially, while inside the well it is a stationary wave, i.e. a 
superposition of two opposite progressing waves. This is also what we have here: 

- the sum over x means that you have to consider the values of x, solutions of \i(x) = 
Aj which have this property. When Aj belongs to [a», 6»] (the suport of the density of 
eigenvalues of the i th matrix), the equation Xi(x) = Aj has no real solution, it has only 
pairs of complex conjugate solutions, which give the stationary wave. The sum of the two 
complex solutions will give rise to some real expression for ijji- n , involving cosinus and sinus 
functions instead of exponentials (cf When Ai is outside [ai,&i], you have to keep 
only the solution which decreases exponentially at infinty. 

From now on, we will consider only the first case, i.e. Aj G [a;,6;]. 

1.5. Kernels 

Remember that we have introduced the orthogonal polynomials in order to integrate 
the joint-density (|1.2|) over a subset I of the N x (p + 1) variables. In this purpose, let us 
as usual || rewrite the Vandermonde determinants: 

A(A ) = T[(Xo;k - A ;i) = det (\ ;i) h 

-*- k,l 

k<l 

Since linear combinations of columns preserve the determinant, we can rewrite: 

A(A ) = cte detVk(Xo-i) 

k,l 

A(A P ) = cte det Vk(X p -i) 

k,l 

The cte is a normalization which comes from the fact that the polynomials Vk and Vk are 
not monic (actually cte = n^o 1 &( n ) N ~ 1 ~ n )- Any partial integration of ( |1.2|) can thus be 
written as an integral over the ijji- n and ifjj. m . Since they are orthogonal, the integration is 
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easily performed, and the final result can be written in terms of 2(p + l) 2 Kernels denned 
by: 

1 N ~ 1 



(1.19) 



and 



Xj) — < 



n=0 



( o if i > j 

-f[2c(X i -X i+1 ) 2 + V i (X^ + Vi +1 (\ i+1 )] if i = i _! 

3-1 3-1 

it 

^ " l=i+l 



, 3-1 3-1 

/ II dXl II ^,m(Ah^+0 



if % < j - 1 



1.20) 



In the p = 1 case discussed in there were only four kernels (the Kij), indeed, the Eij 
which were just numbers were absorbed into the normalizations. But in the general case, 
the contain integrations and cannot be absorbed. Note that the are the propagators 
from to ifjj (i < j): 

dXj i>j;n{Xj)Eij(\i, Xj) = 4>i;n(Xi) 



dXi ipi. n (Xj)E iyj (Xi,Xj) = ipj-n(Xj) 
We thus have the following projection relations: 

dXj EijEji =Eu if i < j <l 



J dXj EijKij =Ku if i < j 



dXiEijKu =Kji if i < j 



(1.21) 



dXjKijKji =jjK u 



1.6. Correlation functions 

In terms of these kernels, the joint density ( [Of ) of all the eigenvalues of all the matrices 
can be rewritten: 

p = cte det Kq >p det -Eo.idet E\^ . . . det E p _i jP 

To obtain the densities and correlation functions of some set of eigenvalues, we have to 
partially integrate p with respect to the other eigenvalues, and this can be done || with 



9 



the help of the projection rules (1.21) . The general result is given in appendix A. Here, 
we will only consider the one and two-points functions. 

The density of eigenvalues (the one-point function) of the z th matrix is: 



Pi(Xi) = K iii (X i ,X i ) (1.22) 

and the two point connected correlation function of one eigenvalue A^ of the i th matrix , 
and one eigenvalue [ij of the j th matrix is: 

p£-( a ^ N) = -Ki,j[Xi, //J/\,.,(//,.A,) + jyK id (\i, fj,j)Eij(Xi, fij) (i < j) (1.23) 

We now have to evaluate the kernels K^j and in the large iV limit. The first step 
will be a generalization of the Darboux-Christoffel theorem, which allows one to rewrite 



Kij as a sum of a small number of terms, instead of the sum of N terms as in (|1.19|) . The 
second step will be to use the WKB approximations (1.18 ) for the ^'s. The kernels Eij 
will be evaluated by steepest descent. 

1.7. Generalization of the Darboux-Christoffel theorem for the kernels 

As in || the Darboux-Christoffel theorem can be generalized. Formally, we write that 

VV>(Ao) = £- (7V - n) ^;iv(Ao) , i>j;n{\) = X N -% ;N (X p ) (1.24) 

and we sum up the geometrical series in Ql.l9| ) : 



Kij(Xi,Xj) = l^_^ ;iV (A l (x))^ ;Ar (A J (y))) (1.25) 

(we have called y the operator x acting on the second variable) . Multiplying both sides of 
( P..25Q by (Xi(x) — Xi(y)) would give in the left hand side a differential polynomial acting 
on (indeed A^(y) can be rewritten as a polynomial in Xj and Pj with the help of 
eq.( |1.7|) ,( |1.9|) ), and in the right hand side a polynomial in x, and y, i.e. a small number of 
ipi- n and ipj;n with \n — N\ N, (an explicit example is given in Appendix C). However 
we shall not do it, but use directly ( 1.25 ) in the large iV limit, where x and y become 



numbers x and y. 

The kernels can thus be approximated by: 



1 y 

Kij(Xi,Xj) ~ — ipi. N {Xi{x)) ^j. N (Xj{y)) 

IS x — y 
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and using the WKB asymptotic expressions (1.18) of i/j^n and ipj- t N : 



k *j ~ E 



/ n x ^ 1 1 1 -iVc (/ x (A l+1 - Ai_i)A- - J y (X J+l - Aj_] 



- 



■<■</ 



ViVcV 2z7riV x - y Jx^xjx^y) 

(1.26) 

where as usually, the x and y stand for the complex solutions of Xi(x) = Aj and Aj(y) = Xj. 
One can also find an asymptotic expression for the kernel E{j by steepest descent: 

where 

0y(A, aO = J] - Az+i) 2 + ^VJ(A0 + ±-V l+1 (X l+1 )) 
i=i ^ ' 

where the A/ (A, /z) are determined by the saddle point equation: 

A* = A , A J=/ u , 2Ai + — Vi(Xi) = Xi-x + X l+1 for % < I < j 



and Dij is the determinant of the matrix of second derivatives of Uij with respect to the 



A,'s: 



= det 



2 + 



V£i(*i+i) 



2 c 



-1 





-1 ... 

2 + ^g^ -1 








-12 + 



Vj- 2 (>j-2) 

2c 



... 

In the particular case x = y we have: 



2 + 



y/-i(^-i) 



2 c 



Uij(x,x) = -J (Aj+i - Ai_i)A- + y (Aj+i - Aj_i)A^ 



and 



-^ij (^7 



^ a:(x)a;.(x) 

t^AJ(x)AJ +1 (x) 



Substituting (|1.26|) and ( |1.27| ) into ( |1.22|) and ( |1.23| ) we can now evaluate the corre- 
lation funtions. 
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1.8. Correlation functions in the short distance limit 



case 1 = 3 



Setting % = j in (|1.26| ) gives: 

1 1 



2t^a;-y > /AJ(x)AS(y) 



where x and y are the complex solutions of Ai(x) = A and \i(y) = \i. When A is close to 
H, at leading order, we keep only the values of x and y such that \x — y\ is small, then this 
reduces to: 



1 1 



sin{iVc(A - /i)Im(A i+ i - Ai_i)} e 



ttN X- 

In particular, when A = \i we obtain the density: 



c 1 

Pi (A) = Kii(X,X) = Im(A i+ i(x) - A;_i(x)) = ImPi(x) 

TV TV 



(1.28) 



When A is close to \i but different, we can compute the two-point connected correlation 
function: 

pf}{^ y) = -Ki,i(X,fi)K iti (n, A) 



i.e. 



( sin Ntt(X — n)pi(X) 

A— V A^7T(A - jl) 



(1.29) 

we recover the universal two point correlation function in the short distance regime. 
• case i ^ j 

It is now meaningless to consider the limit A^ close to Xj since they are eigenvalues 
of different matrices. Generically, is of order 1/N, which means that the connected 
correlation function is of order 1/N 2 , and we can say that in the large iV limit Aj and Aj 
are uncorr elated. 

The only limit in which the correlation may become larger than 1/N 2 is the case 
where x — y becomes small. The equation x = y defines a function Xi(Xj). The problem 
is that this function takes complex values in the interesting domain - for example we see 
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from eg ( |Q||) that Aj+i(Aj) or X i _i(X i ) take complex values, ( this fact has already been 
debated in || ) - and we have not found any physical interpretation, except in the case of 
the continuous model described in the next section. 

However, let us assume that \x — y\ is small (we will also assume i < j). We introduce 
the scaling variable: 

A 



In the limit \x — y\ — > 0, the Taylor expansion of the term appearing in the exponential in 
( |L26|) gives: 



ij 



(A i+ i - A,-!)^ - / (A i+ i - Xj-i)^ = -Uij + '-jj^(x-y) 2 + 0((x - y) 3 ) 



Therefore, we have: 



/ 7T \ V 1 rc~ +NcU l3 1 -A 2 
Kij ~ ( — I != 1 1 e — e 



ViVc/ 2inVN\ Dy A 

Kj , . ( - ) - 1 ^e~ NdJ " Ie A2 (1.30) 
3 \Nc) 2inVN\ Dij A 1 ' 



1 /m 2 ? 1 rc~ -NcU,, 



—E tj ~ [ — ] 



AT tJ \NcJ ^VN\ 
Remember that each of these expressions is actually a sum over the different values of x 
and y satisfying \i(x) = A and Xj(y) = fi- But since A and [i are not assumed real, the 
different values of x or y contributing to this sum are not the complex conjugates of each 
other, and the exponential term cannot simply be replaced by a sinus function. 
However, we have the correlation function: 




(1.31) 

To go further one would have to make some asumption on the argument of x and y (i.e. 
where in the complex plane are A and n), and then separate the imaginary and real parts 
of Uij and A. We would then observe that in this scaling regime (\x — y\ ~ 1/y/N), pfj 
is of order 0(1/N) instead of 0(l/iV 2 ), and that it presents some kind of universality: it 
doesn't depend very explicitely on the potentials V^(A). The signification of this correlation 
function is therefore not clear, and has to be further studied. However, this calculation 
was performed to prepare the next section, the continuous chain of matrices, where it is 
possible to have both A and \x real and (x — y) small, because will be small too. 
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1.9. Smoothed functions 

If we smooth expression ( |1.23| ), (we keep in the sums over x and y only the terms 
which don't oscillate at frequency N, i.e. for which the exponentials exactly cancel), we 
obtain: 

Plj(Xi, Xjhmooth ~ - 47r2JV 2 (x(A^) - y(Xj)) 2 com P lexcon J u S ate ( L32 ) 
This result is exactly the same as for two matrices |10|,|S| . 



2. continuum limit 

The model of a chain of matrices is naturally extended to the model of a time depen- 
dent random matrix. We replace the integer index i by a continuous time t = ie, wich runs 
from to T = pe. The coupling term c(Mj — Mi + i) 2 becomes a kinetic term, we set: 



1 




Vi(X) — ► eV(X,t) 

The partition function thus becomes a functionnal integral for a one dimensional matrix 
field: 

f -Ntrtfdt (i(f) 2 + W),i)) 

Z= D[M(t)] e V J 

We will now just make use of the results of the previous section for the discrete time 
case, and translate them into the language of the continuous model. 

We now have a set of orthogonal wave functions depending on time 

tpi-n{X) — >ip n (X,t) and ^i-n(X) — > ^ n (X,t) 
which satisfy the orthogonality relationship: 

J &Xll) n {X,t)$ m {X,t) = S n;m 

Actually, one has to change the normalization, just to absorb the constant factor (tc/Nc) 1 ^ 2 
of eq (1.18) which becomes infinite in the continuum limit. But this doesn't change any- 
thing else. 
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2.1. The function X(x, t) 

We can define a family of functions X(x,t) of an auxilary variable x: 

Xi(x) — ► X(x, t) = J2 (*k(t)x k (2.1) 

k 

Note that acording to eq ( P-.13|) , k runs from — oo to +00 (except for the gaussian case), 
so that ( |2.1p has to be taken as a formal expansion in Laurent Series. 
X(x,t) obeys the continuous limit of the equations of motion ( |1.15|) : 

d_ 2 

at 

with the boundary conitions ( |1.16| ): 



2 X(x,t) = V'(X(x,t)) (2.2) 



X(x,0) =ax + ^2a k (0)x~ k ; X(x,0) = — — + 0(l/x 2 ) 

k<0 

X(x,T) = -a + y"x k a k (T) ; X(x, T) = — + Oix 2 ) 
x f— ' a 

k>0 



those last conditions can be rewritten in a way which doesn't involve any expansion in 
powers of x ***: 

X(x, 0) ~ — — - — — when A(x,0)^oo 

X(x,T) ~ — — — when X(x,T) — > 00 
X(x,T) 

however, to define the function x(X, t) you can't avoid performing the formal expansion. 



Remark: eq (|1.14|) becomes: ka>kOt-k = 1 



2.2. The momenta 

The momentum P(x,t) is the time derivative of X(x, t) (eq ( |1.7| ) and (|1.8|)): 

dX 



p(i) = A at 



- v(X, t) — z7rp(A, t) 

its imaginary part p(X, t) = — ImP(t) is the density of eigenvalue at time t. 



*** Actually, even in the discret case, cjo(Ao(^)) = 2c(Xo{x) — Ai(x)) + Vo(Ao(a;)) is the resolvant 
of the first matrix u>o(X) = jj < Tr A _ 1 M >, and oj p (X p (x)) = 2c(X p (x) — X p -i(x)) + V p {X p {x)) 
is the resolvant of the p th matrix. In the continuum limit we have lu(X(x, 0), 0) = — X(x, 0) and 
uj(X(x,T),T) = X(x,T). For the intermediate values of i (intermediate times), there is no simple 
relationship between the resolvant cJi(A) (u>{X,t)) and the functions Xi(x) (X(x,t)). 
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2.3. The Kernels 

The discrete kernels K^j and E^ naturally have some continuous version: 



N-l 



K it j(X, fi) — > K(X, tin, t') = — ^2 ^n(A, *)$ n (/*, t' 



n=0 



A(f)=/x -Nj!'dT (lA 2 + l/(A(r),r; 

D[X(r)} e V 

A(i)=A 



Eij(X,fi) E(X,t\fjL,l/) 
E is the usual quantum mechanics propagator for a single particle in the potential V: 



1p n (X,t') = j dpip n (p,t) E(/J,,t\X,t') 

And those Kernels allow to calculate the two-point connected correlation function (t < t')\ 

pV(\,t\n,lf) = -K(X,t\fi,t') (K{(iJ\\,t) - ^E(X,t\p,t'yj 
2.4- Correlations 

Let us calculate the 2-point correlation in the limit where 5X = X — p and St = t — t' are 
small, of order 1/N. In this purpose, we have to evaluate the kernel K(X, t\X+8X, t+5t) with 
the WKB approximation (|1.26| ) (remember that we have absorbed the (tt/Nc) factors): 



K ~ 2Re 



1 1 _ x -N (j X P(X',t)dX' - J X+SX P(X',t + 5t)dX' 



2mN sfxTX 1 8x 



where -P(A, t) = A = (dX/dt) \ x = v — iivp. Writting SX = X'Sx + PSt the denominator 
X'Sx is 5X — vSt + mpbt, and the term in the exponential is P5X + 5t j X (dP/dt) \\ which 
is evaluated through a play with the partial derivatives at A fixed or x fixed: 



V'(X) = X 



dP 
~dt 



dP 
~dt 



+ 



dP 

ax 



dX 

In 



We thus obtain: 
K ~ Im 



1 



Nn(v5t -8X- iirpSt) 
Similarly, observing that the continuous limit of Dij is 

-Da ~ 28t 

c 



Niirp(v8t-8X) N(v5X + 5t(V - \v 2 + \n 2 p 2 )) 
e e (2.4) 
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the asymptotic expression (|1.27 ) for E turns into: 



l E 1 + V8t) 



for St > 



N y/2irN5t 

• equal time correlations: 

When St = we recover the usual one-matrix-model result: 

T ^ sin (NnpSX) NvSX 

K ~ e 

NnSX 

For instance when SX = we obtain the one-point function (the density) K 
when SX ^ but of order 1/N, we have the universal correlation function: 



(2.5) 



p. And 



> (c) (A,t | X + SX,t) 



-P 



sin (NirpSX) 
NwpSX 



So, at equal time, the correlation function of close eigenvalues is still universal. 
• different times: 

When St 7^ 0, the denominator of K never vanishes. It reaches its minimum for 
SX = vSt. We shall thus consider the regime: 



where we have: 



and 



K 



SX = vSt and St ~ 1/N 



p NSt(V+±v 2 + ±n 2 p 2 ] 



NTT 2 p 2 St 



:E 



p^/n -NSt(V + \v 2 ) 



N ^/2NK 2 p 2 St 



Thus 



1 I 1 /7T e 



p (c) (A,t | X + vSt,t + 6t) ~ -p 2 — \ — + 



A A V 2 VA 



|A\ 



with A = Nn z p z St 



(2.6) 

which is remarkably universal, it doesn't depend on the precise form of the potential 
V(X,t). 
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3. Conclusion 

We have thus, as announced in ||, generalized the properties of the correlation func- 
tions of the two matrix model to the chain of matrices. 

As expected, the chain of matrices obeys the same universal properties as the one and 
two-matrix models. We obtain similar results: 

- the correlation function is generically of order 1/iV 2 , except in the regime where the 
denominator of the kernel K becomes small. This may happen only at small distance SX 
and at small time interval St. 

- At equal time (St = 0) we recover exactly the well known universal correlation function 
of the one-matrix model. This is therefore one more matrix model where this universality 
is proved. 

- At different times (St > 0), the 2-point connected correlation function is again a universal 
function of St and the shifted interval SX — vSt, i.e. it doesn't depend on the potential 
V(X,t) (indeed the term involving V exactly cancels when multiplying (|2.4| ) and (|2.5| )). 

Moreover, we have given a very explicit way to compute the correlation functions. 
This method (the generalized Darboux Christoffel theorem) might allow one to compute 
the correlation functions even for N finite, provided one knows the orthogonal polynomials. 
What we have computed here is an approximation in the large iV limit, using the WKB 
approximations of the orthogonal polynomials. This method has shown that the function 
x(X) plays a very important role, as the best auxilary variable, its physical meaning should 
be better understood, beyond the orthogonal polynomial's method. 

For instance, another way of studying the chain of matrices is through the loop equa- 
tions (the invariance of the partition function w.r.t the change of variables, see Appendix 
D), which would allow to compute the following subleading terms in the 1/N expansion. 
The loop equations in the large iV limit give some algebraic equation (of very high degree) 
for the resolvent of the first (or last) matrix of the chain. The parametrization of A and 
u(X) as functions of x allow to solve this algebraic equation. 

Acknowledgements 

This work was partly supported by the TMR Network contract ERBFMRXCT 960012. 



18 



Appendix A : Correlation functions of any set of eigenvalues 

The general correlation function of any set of eigenvalues has been computed by || . 
In addition to the kernels and Eij, we will need to introduce the H^s : 

Hij = Kij - jr E ji if 3 < i > H ij = K ij if 3 > i 
To compute a partial density of eigenvalues, for instance 

Pn ,...,n p (A();l, ■ ■ ■ i ^0;n i • • • ? A p; i, • • • ? ^p;n p ) 

you have to consider all the possible permutations of all the n = uq + . . . + n p variables, 
and then (with obvious notations): 

p(Ai, . . . , A„) = det H iyj = ^(-l)' 7 Yl H ii<T ^(Xi, A^)) 

a 1 

The connected correlations are obtained by reducing the sum to the cyclic permutations 
only. 

For instance, we have (we assume i < j < k): 

Pi(Xi) = Kij(\i, A») 
Pij (Ai, Xj) = —Ki t j(Xi, Xj)Kj i(Xj, Aj) + —Ki t j(Xi, Xj)Ei t j(Xi, Xj) 
Pij,k = KijKjkKki+Ki k K k jKji- — (KijKj k E ik +Kj i K ik Ej k +K ik K k jE i j^ 

Appendix B : Explicit expressions of orthogonal polynomials 

The orthogonal polynomial V n and V n have explicit expressions as matrix integrals: 

-p m (am am M\T\ ~ NtrV i( M i) tt ~N ctr (Mi - M i+1 

V n [X) = / dM 0nXn ...dM pnXn det(A-Mo) [[ e [[ e 

i=0 i=0 

V n (X)= dM 0nXn ...dM pnXn det (X-M p ) l[ e lK l> J] e V 

2 = 1=0 
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Appendix C: An explicit example: the gaussian potential, V(A) = \g\ 2 



3.1. Discrete case 

Let us consider a constant gaussian potential: V^'(A) = gX. We will use the following 
parametrization g = 4c(coshx — 1) for later convenience. 

The operators X\ contain only two terms, which we write as: 

c „ [ft n [ft 1 

The equation of motions are linear and thus easily solved: 

a i+ i + ati-i - 2ai = ^-cti , (3 i+1 + - 2/3^ = ^-/^ 

2c 2c 

with the boundary conditions: 

a = /3 p = a , a 1 = (l+^-)a , j3 p -i = (1 + ^-) a 

{1 + Yc )ap - ap - 1 = 2^ 
The first line implies that $ = a v -i. The solution is then 

cosh(z + i)x 2 1 

«i = a t-^ — where a = — — t- 

cosh 4c smh (p + l)x tanh ^ 

The momentum Pi = is: 

. 1 /n 1 A „1 + An. 1 ^ In 





I 1 - -2-^Aj \ \/-- , P l = -2—^X t -} - -.v 

where we have introduced Ai = sinh (p — 2i)x/ sinh (p + l)x and 

cii = 2 v /a i a p _i 

which will play a very important role as the natural scale for Aj. Actually 2a^ is the width 
of the distribution of eigenvalues of the i th matrix. Indeed, in the large N limit and n ~ N, 
we can eliminate x and write Pi as a function of A^: 



2 . Xi A? 



Pi{\) = -- [A^-.j^-l 
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its imaginary part is the density of eigenvalues of the i th matrix: 



P*(A*) = ^/l-^ 



which is a semi-circle law of diametre 2aj. 

The wave functions can be exactly computed in terms of Hermite polynomials H n (x): 



1pi;n{Xi) = Ci [ -J 



a n -i \ 2 1 



, -X\ -Nil-Ai) 
H n e 
n' V a 



ft 



H n ( 2^^- I * 

a. 



(where Cj and Cj are some constants, irrelevant for what we need, we just know that 

1 / 2JV \ 



Ci Ci 



«i V T 



And, we can write some WKB approximations in the large A limit 
VAc/ sin ^ \ ctj / 





2n-2iV + l 




\ 4 




COS 







1 NJ A 7T 

(n+ -)0i - y sin 20; - - 



AA 



A? 



~ / 7T \ f"# 2(2tt)^ / aj 

^i;»(Ai) ~ TT / ■ , 

Vive/ A/ajSin^j \otp-j 





4 




( ) 




COS 


\a p -jj 







(n+-)(f) j - — sm20 j - - 



-AA ? -4 



where A, = a, 



cos ■ 



Now let us compute the kernels: 
The propagator Eij is: 



7 — i — 1 

/ 7T 



sinhx 



~ Nc ^J^ [^sh (j - i)x + AJ) - 2A,A,] 



sinh (j -i)x 



This is an exact result since the saddle point method is exact in the gaussian case. We 
also have the determinant of the second derivatives of the potential Dij(Xi, Xj) which is a 
constant: 

= sinh (j-i)x 
tJ sinh x 

The Kernel Kij(Xi, Xj) obeys a generalization of the Darboux Christoffel Theorem: 

A,- _ A 



2 ^ _ 2r^ + ( r - -) (1 + A/) — + 
a,- a, r V a, 



A, a, 9 



a,- 2A dXj 
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/ eg 

C^p—i 



'j;N-V 



p—i 

OLi 



at dj r V a i 2N d\i J ) ]/ a p-j V a i 



where r denotes the ratio r = \Jct p -iCtj/cti(y. v -j. Let us emphasize that these equations 
are exact, even for iV finite (remark that when % = j we have r = 1 and we recover the 
usual Darboux-Christofel theorem). With the operatorial notations, they are obvious, they 
are just the rewritting of 

y 

{Xi(x) - Xi(y)} Kij = {Xi(x) - Xi(y)} ^—r^i-N^j-N 

x y 

{Xj(x) - Xjiy^Kij = {Xj(x) - Xj(y)} -t^^n^n 

x y 

(if you want to rederive them, be carrefull that here, the partial derivative with respect to 
Xi carries on the bra, it is aP-). 
In the large N limit, we have 



/ 7T \ V 1 N(Ai cos 2 - Aj cos' 

~ I — I — - e 

V Nc J yj en sin <pi aj sin <pj AnN 



(1 + J) sin \{(Pi - (Pj) siniV(Ci - 0) + (I-7) cos - fa) cosN(Q - Q) 



1+T S 

2r 



COS {(pi - <pj) 



+ (same thing with <pj — > —(pj) 
where we have noted ((fa = (p — \ sin 2<p. 
- For % = j we have r = 1 and we recover: 

/ \ \ r Ai(\2 ,.2\ 

1 1 /sinAr(C(A)-C(/x))\ "2" of ' ^ 



4yriV 



sin 



which gives for A = the density of eigenvalues: Pi(X) = sin</>. 

- When i 7^ j, we remark that the denominator never vanishes. It is maximum when 
(pi = (pj, i.e. Xi/cii = Xj/cij, this is the regime in which one could have a not too small 
correlation. 
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3.2. Continuous limit 



It can be obtained by setting c = ^, t = it, T = pe, g = ez/ 2 , and taking the limit 
0. 

We find: 



a 2 



v sinh vT 



which we write 



a(t) is given by 



1 

A(x, t) = a(x cosh vt H — cosh v\T — t)) 

x 



A(0, i) = a(t) cos(p 



a(t) = 2aA/cosh vt cosh u(T — t) 



In the small distance regime (5X small and 5t small), we have (up to some factors 
which will cancel when we compute the correlation function): 

11 5(j) sin {Nti pa sin 050) + -^St cos (Nirpa sin 

K(X,t\X + 5X,t + St) - 



2irNasm(j) 2ff + \S(j) 2 



b 2 

. ~FT IW dill U/ r, i — 

—E(\,t\\ + 5\,t + 5t)~ - e 2 

y/2ivN6t 



1 , 1 -f ( a 2 sin 2 0^ + * cos 2 #t) 



And thus we have the two-point connected correlation function: 

p (c) (A, t\X+SX, t+St) = -if (A, t|A+5A, (a:(A + SX, t + St\X, t) - ^E(X,t\X + SX,t + St) 
which gives: 

^ 1 5(f) 2 sin 2 (Nn pa sin <j>5<j)) — -^St 2 cos 2 (Nnpa sin </><50) 

P 4n 2 N 2 a 2 sin 2 (2^ + ±50 2 ) 2 

^ 1 <5t 50sin(A/>pasin0<50) + ^<5tcos(A%pasin0<50) ^ ^ 



(2yriV(St) 3 / 2 asin0 2^ + ±50 2 

■f (a 2 sin 2 0^ + ^ cos 2 05t) 



c 



Appendix D: The loop equations 
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The loop equations are just the consequence of the invariance of the partition function 

-iVtr Vi(Mi) + c £,(M, - M l+1 ) 2 



J dM 



.dM p e 



with respect to some change of variables Mi — ► Mi + e/(M , . . . , M p ). They are most 
simply written by introducing the resolvant 



N \ z-M 

Actually, one can easily obtain close equations in the large N limit only for the resolvant 
of the first (i = 0) and last (i = p) matrices. We will give here the recipe to write 
the algebraic equation for ui{z) but without any proof (the proof obtained by multiple 
recurrences is quite tedious but not complicated): 
first we define the functions: 

1 1 

z (z) = z , Zl (z) = z + 2^Vo(*0 - Yc U ^ 

and z i+1 (z) = 2zi(z) + —V- (zi(z)) - Zi-^z) for 0<i<p 

each function Zi(z) is a polynomial in z and ui{z). The algebraic equation relating z to 
w(z) is then: 

(z (z) - Zl (z) + ^-V^zo(z))^ (z p (z) - zp-^z) + y c V p(z p (z))^ = ^(*o(*), • • • > 

(3.2) 

where the right hand side P(xq, . . . polynomial of p + 1 variables (now considered 

as independent): 

P(x , ...,x p ) = Pol X0 _ Xp 1 (tr f p+l (x , . . . , x p )^-^ • • • 

and the functions f n (xo, X2, ■ ■ ■ , x n -\) are defined as follows: 

/o = l , fi(x) = x + ±V£(x) 

fn+l(X(), ■ ■ ■ , X n ) = f n (xo, . . . , X n -i)(x n + — V^(x n )) — / n _i(x , . . . , X n -2)X n X n -i 

Note that all this is symetric in the exchange of the two extremities of the chain i — > p — i. 
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Most of the coefficients of the polynomial P are unknown, they are some arbitrary 
constants which have to be fixed by some analytical considerations about ui{z). 

- One of the constraints is that u(z) ~ 1/z when z — > oo. 

- The other constraints require more physical input. For a single well potential Vi(\), 
we expect the density of eigenvalue to have a connected support, i.e. that u(z) has only 
one-cut. 

This one-cut asumption should allow in principle to determine all the unknown coef- 
ficients of P. We made exactly the same asumption when we replaced the operator a; by a 
number x in the large iV limit, and we observe that the functions Xi(x) give a parametriza- 
tion of the solution of eq ( |3.2|) . Actually, if z = Xo(x) we have that each Zi(z) is nothing 
but Xi(x), and the loop equation ( |3.2|) is nothing but the product of (|1.16|) and fll.lTf ). 
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